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Abstract 

Let u = {u(t,x),t € [0, T],x £ M rf } be the solution to the linear stochastic heat 
equation driven by a fractional noise in time with correlated spatial structure. We study 
various path properties of the process u with respect to the time and space variable, 
respectively. In particular, we derive their exact uniform and local moduli of continuity 
and Chung-type laws of the iterated logarithm. 
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1 Introduction 

Stochastic analysis of fractional Brownian motion (fBm) naturally led to the study of 
stochastic partial differential equations (SPDEs) driven by this Gaussian process. The 
motivation comes from wide applications of fBm. We refer, among others, to 
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|3U] and [33] for theoretical studies of SPDEs driven by fBm. To list only a few examples 
of applications of fractional noises in various areas, we mention [15] for biophysics, [5] for 
financial time series, m for electrical engineering, and [7] for physics. 

The purpose of our present paper is to study fine properties of the solution to the 
stochastic heat equation driven by a Gaussian noise which is fractional in time and colored 
in space. Our work continues, in part, the line of research which concerns SPDEs driven by 
the fBm but at the same time it follows the research line initiated by Dalang [8] which treats 
equations with white noise in time and non trivial correlation in space. More precisely, we 
consider a stochastic linear equation driven by a Gaussian noise that behaves as a fractional 
Brownian motion with respect to its time variable and has a correlated spatial structure. 
A necessary and sufficient condition for the existence of the solution has been given in [3j 
and other results has been given in mi m m m among others. 

To briefly describe the context, consider the stochastic heat equation with additive 

noise 

c) 1 } I 1 

— = -Au + W, te[0,T], xeR d , (1) 

u(0, x) = 0, x £ M d , 

where T > 0 is a constant and W is usually referred to as a space-time white noise. It 
is well-known (see for example the now classical paper [8|) that dU) admits a unique mild 
solution if and only if d = 1. This mild solution is defined as 

u(t,x) = f f G(t — s,x — y)W(ds,dy), f € [0,T], iGK. (2) 

Jo J R d 

In the above, W = {W(t, A)-,t > 0, A € £>t,(M d )} is a zero-mean Gaussian process with 
covariance given by 

E (W(t, A)W(s,B)) = (s A t)X(A n B), (3) 

where A denotes the Lebesque measure and Bb(^ d ) is the collection of all bounded Borel 
subsets of M. d . The integral in ([2]) is a Wiener integral with respect to the Gaussian process 
W and G is the Green kernel of the heat equation given by 


G(t, x) 


(2irt)~ d / 2 exp if t > 0, x G R d , 

0 if t < 0, x G M 01 . 


(4) 


Consequently the mild solution {u(t,x),t e [0,T],x € M} is a centered two-parameter 
Gaussian process (also called a Gaussian random field). It is well-known (see e.g. [33] or 
mi that the solution defined by the ([2]) is well-defined if and only if d = 1 and in this case 
the covariance of the solution ([2]), when x € M is fixed, satisfies 

E (u(t, x)u(s, x)) = ^ y/t + s — y/\t — s \j , for every s, t € [0, T\. (5) 
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This establishes an interesting connection between the law of the solution m and the 
so-called bifractional Brownian motion introduced by M- Recall that, given constants 
H G (0,1) and K G (0,1], the bifractional Brownian motion (Bf ,K )te[o,T] is a centered 
Gaussian process with covariance 

R H,K {t,s) := R(t,s) = ({t 2H + s 2H ) K - \t-s\ 2HK ^j , s,le[0,T]. (6) 

Relation ([5]) implies that, when x G R is fixed, the Gaussian process {u(t,x),t G [0,T]} is 
a bifractional Brownian motion with parameters H = K = -( multiplied by the constant 
2~ k —)=. Therefore, many sample path properties of the solution to the heat equation 

v 27 T 

driven by space time white-noise follow from mmmm- 

In order to avoid the restriction to the dimension d = 1, several authors considered 
other types of noises, usually more regular than the time-space white noise. An approach is 
to consider the white-colored noise , meaning a Gaussian process W = {W(t,A)-,t > 0, A G 
£>;,(M rf )} with covariance with zero mean and covariance 


E(W{t,A)W(s,B)) = (t As) [ [ f(z-z')dzdz'. (7) 

J A Jb 

Here the kernel / is the Fourier transform of a tempered non-negative measure /x on M rf , 
i.e. /x is a non-negative measure which satisfies: 


1 


/x(d£) < oo for some l > 0. 


.i + iei 2 . 

Recall that the Fourier transform / of /x is defined through 


[ f(x)ip(x)dx= I r<p(SMdt), VvJGS(R d ), 

JR d JR d 


( 8 ) 


where 5(M d ) denotes the Schwarz space on M rf . A useful connection between the kernel / 
and its associated measure /x is the Parseval inequality: for any if, vJj G 5(M d ), 

/ / <p(x)f(x ~ y)i>{y)dxdy = (2vr) -d f ^(^^(Om^O- ( 9 ) 

JR d JR d JR d 

In this work, we will include the following two following basic examples. 


Example 1 The Riesz kernel of order a: 

f(x) = R a (x) := 7 a,d\x\~ d+a , (10) 

where 0 < a < d, = 2 d ~ a Tr d / 2 T((d — a)/2)/T(a/2). In this case, ^(df) = 
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Example 2 The Bessel kernel of order a: 


poo 

f{x) = B a {x) := y a / u ,(«-‘0/2-i e -« e -|*l n^w)dw, 

Jo 

where a > 0, r y' a = (47r)“/ 2 r(a;/2). In this case, /r(d£) = (1 + |£| 2 ) _ “/ 2 d£. 


For any function g E L 1 (R d ) we denote by Tg the Fourier transform of g, i.e., 
(Bgm= [ e- i ^g(x)dx, £ € R d , 

JR d 

where (•, •} denotes the inner product in R d . Recall that for any t € R and x € R d , the 
Fourier transform of the fundamental solution Q is given by 

FG{t,x - •)(£) = exp !{*>()} (£)> CeR d . (11) 


Dalang [ 8 ] proved that the stochastic heat equation with white-colored noise admits 
a unique solution if and only if 

[ < 00 • ( 12 ) 
Jr* i + |£r 

Obviously, this condition allows x to be in higher dimensional space. For example in the 
case of the Riesz kernel, the stochastic heat equation with white-colored noise admits an 
unique solution if and only if d < 2 + a. 

Under (j!2j) . the solution of CD) with white-colored noise can still be written as in 
(| 2 |). One can compute the covariance of the solution with respect to the time variable (see 
e.g. [28]). For fixed x E R rf , and for every s < t it follows from ParsevaPs identity and ED 

E (u(t,x)u(s,x)) = {2n)- d f du f g,(d£)e~ 2 (*— u )l^l 2 e — 2 ( s_ , ( 13 ) 

Jo Js. d 


In the case when / is the Riesz kernel (i.e., g(df) = |£| a df, we get 

E(u(t, x)u(s, x)) = (27 T)~ d jf ^ (it + s) 1 - V - (t - s) 1 'V) . 

In this case, the solution coincides, modulo a constant, with a bifractional Brownian motion 
with parameters H = ^ and K = 1 — Thus the sample path properties of this process 
can be deduced from nsi eh es]- 

Our purpose in this paper is to study the linear heat equation driven by a fractional- 
colored Gaussian noise (see section 2 for details). When the structure of the noise with 
respect to the time variable changes and the white noise is replaced by a fractional noise, 
the solution does not coincide with a bifractional Brownian motion anymore (see [ 6 ] or 
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[55]). Some new methods are needed for analyzing the path properties of the solution 
with respect to the time and to the space variables. By appealing to general methods for 
Gaussian processes and fields (cf. e.g., [2D] SS r 4D]), we establish the exact uniform and 
local moduli of continuity and Chung type laws of the iterated logarithm for the solution 
of ([T|) with a fractional-colored Gaussian noise. 

Throughout this paper, we will use C to denote unspecified positive finite con¬ 
stants which may be different in each appearance. More specific constants are numbered as 

Cl, C2,.... 

2 The solution to the stochastic heat equation with fractional- 
colored noise 

We consider a Gaussian held {W H (t,A),t € [0,T],A € £>(,(R d )} with covariance 

E(W H (t,A)W H (s,B)) = R H (t,s) [ [ f(z-z')dzdz', (14) 

J A JB 

where Rn{t,s) := \{t 2H + s 2H — \t — s\ 2H ) is the covariance of the fractional Brownian 
motion with index H and / is the Fourier transform of a tempered measure /r. This noise 
is usually called fractional-colored noise. We will assume throughout the paper that the 
Hurst parameter H is contained in the interval (^, 1). 

Consider the linear stochastic heat equation 

I 1 

— = -Ail + W h , te[0,T],xeR d (15) 

with vanishing initial condition, where {W H (t, x ), t £ [0, T] , x € is a centered Gaussian 
noise with covariance (1141) . In the following we collect some known facts: 

• A necessary and sufficient condition for the existence of the mild solution to the 
fractional-colored heat equation (1151) has been given in [3] (see also ESI). Namely, 
(USD has a solution {u(t,x),t >0,i£ R d } that satisfies 

sup E (u(f, x) 2 ) < Too 
te[o,T],xeM d 


if and only if 


1 


2 H 


Ll(d£) < OO. 


d + iei 2 . 

When (1161) holds, the solution to (|15D can be written in the mild form as 


(16) 


u(t, x) = / / G{t — u,x — y)W H (du,dy), f € [0,T], x£ 

Jo J R d 


(17) 
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It follows from [3] or |28] that, for x € fixed, the covariance function of the Gaussian 
process {u(t,x),t € [0,T]} can be written as 


E {u(t, x)u(s, x)) 


oih{ 2 tt ) 


-d 



dudv 


o Jo \ u ~ v 


2-2 H 


j 

J 


li{d£)e 


2 e 


( s - v )\£\ 2 

2 


(18) 


with an = H(2H — 1). In the particular case where the spatial covariance is given 
by the Riesz kernel, the process t t —> u(t,x) is self-similar with parameter H — ^ = r L . 
However, for H € (|, 1), {u(t,x),t > 0} is no longer a bifractional Brownian motion. 

For 0 < a < d, the notation 


~ i er a de, (i9) 

means that for every non-negative function h there exist two positive and finite constants 
C and C', which may depend on h, such that 

c' [ h(m~ a dC< f h(£)ii(d£) <C [ hm\- a dd. ( 20 ) 

J R d JR d JR d 

It has been proven in [33] that, under condition (1191) . there exist two strictly positive 

constants C 2 such that for all t, s € [ 0 , 1 ] and for all x € M rf , 

ci\t — s\ 2H 2 ~ < E (|u(i, x) — u(s, x)| 2 ) < C 2 \t — s\ 2H ~2 . (21) 

Remark 1 The Riesz kernel obviously satisfies m>- The Bessel kernel satisfies m> and 

the constants in ( 1201 ) are C = 1 and C' > 0 depending on h. 

Under Condition m, the condition m is equivalent to 

d < AH + a. (22) 

In the sequel (e.g., Theorem |T] below), we will also use the notation fix ) x g{x) 
which means that there exists two positive and finite constants C and C' such that C' < 

f(x)/g(x ) < C for all x in the domain of / and g. In the rest of the paper, we will assume 

that (fT9l) is satisfied. 

2.1 Sharp regularity in time 

For any fixed x € M rf , to further study the regularity properties of the sample function 
t i->- u(t,x), we decompose the solution {u(t,x),t > 0 } in (fT71) into the sum of a Gaussian 
process U = {U(t),t > 0}, which has stationary increments, and another Gaussian process 
whose sample functions are continuously differentiable on (0, oo). In particular, when the 
spatial covariance of noise W H is given by the Riesz kernel (I10p . then U = {U(t),t > 0} 
becomes fractional Brownian motion with Hurst parameter 7 = H — We will apply 
this decomposition to obtain the regularity in time of the solution {uft,x),t > 0 }. 
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Motivated by [21] (see also [23 El] for related results), we introduce the pinned 
string process in time {U(t),t > 0} defined by 


U(t ) = 


(G(t — u,x — y) — G(—u, x — y))W H (du, dy ) 


+ / / G(t — u,x — u)W H (du,du). 

I o Jm d 



Note that 17(0) =0 and U(t) can be expressed as 

u (t)=[ [ (G{(t - u)+,x-y) - G((-u)+,x - y))W H (du,dy). 


(23) 


In the above, a+ = max(a, 0). The following theorem shows that {17(f), f > 0} has station¬ 
ary increments and identifies its spectral measure. This information is useful for appleying 
the methods in [20, 39] SD] to study sample path properties of {U(t),t > 0}. 

Theorem 1 The Gaussian process {17(f), f > 0} given by A23\) has stationary increments 
and its spectral density is given by 


2(2?r) d a H f /*(<%) 

Mr) - 17 pr 5 - l 


(24) 


Moreover, under condition m, we have 

fu{r) > 


|2tf-4=H+l ' 


(25) 


Proof : For every 0 < s < t, by Parseval’s identity Q and relation dill) , we can write 
E [(U(t) - U(s )) 2 ] 

= E ( f f ( G((t - u)+,x -y)~ G((s - u)+,x - y)) W H (du, dy) 

\J R JR d / 

= (27 r)- d a H f y(dO f [ ( e -§Ct-^)|g| 3 f _ e -f(*-«)|g| a i ) 

JR d Jr Jr \u-v\- 2H \ V 




Let ip(u) = e 2 7 u )l£l“l| t>u i — e u )l£l 2 lr s>u \. Then its Fourier transform is 

F^r) = (M ,T - e ~ i>T )— 


(26) 
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By using again Parseval’s relation Q in dimension d = 1 and (1261) . we get 


E [(U(f) - C/(s)) 2 ] = (2vr )~ d a H / 


dr 


= ( 2 vr) d a H [ fJ,(d£) I 

Jr* J 

= 2(2n)~ d an [ [l — cos((f — s)r)] ■ 

Jr 


— \T(p(t )\ 2 

dr 2 [1 — cos((f — s)r)] 


r 2 , jili 

dr f 


\2H—1 


H(d£) 
t 2 + ^ 


where the last step follows from Fubini’s theorem and the convergence of the last integral 
in ]i(d£) is guaranteed by relation (fTBl) . It follows that the Gaussian process {17(f), f > 0} 
has stationary increments and its spectral density is given by (1241) . 

Under condition (fl9l) . we have 


fu(r) x 



\2H-4=2+1 ' 


This finishes the proof of (1251) and thus the conclusion of Theorem [T] is obtained. ■ 

When the spatial covariance of noise W H is given by the Riesz kernel (1101) . we have 

= a H 2 2H - 1 T(H - r df 

(2vr) d+ ^T(l - H) \t\ 2H ~ 1 Jr* |£|« (t 2 + ^ 

— i) I" dr] 

(2vr) d+ ^T(l - H) \t\ 2H -^ +1 Jr* \rj\ a (l + Jlf) ' 


Therefore, in the Riesz kernel case, {17(f), f > 0} is, up to a constant, a fractional Brownian 
motion R 7 with Hurst parameter 7 = H — . 

Recall that the spectral density of the fBm 7 > 7 with Hurst index 7 E (0,1) is given 
by / 7 (A) = c 7 |A|~( 1+27 ) with c 7 = bin l 7 r 7 T( 1 + 2 7 ) ^ Hence, we have the following corollary. 


Corollary 1 Let U = {17(f), t > 0} be the Gaussian process defined by [23 1) such that the 
spatial covariance of W H is given by the Riesz kernel GW- Then U coincides in distribution 
with CqB 1 with 7 = H — and 


{27r)- d+1 2a H 2 2H - 1 r(H - \) r dr] 

sin (vr(d - S=sl)) T(1 + 2i7 - ^)T(1 - H) Jr* \ v \« (1 + Ml) ' 


Now for every f > 0 we have the following decomposition 


u(t, x) = 17(f) - Y (f) 

















where 

Y(t) = f 
J —O O 

The following theorem shows that the sample function of {Y(t),t > 0} is smooth, 
which is useful for studying the regularity properties of the solution process {u(t,x),t > 0 } 
in the time variable. 

Theorem 2 Let x € be fixed and let [ a,b] C (0,oo). Then for any k > 1 there is a 
modification of {Y(t),t > 0} such that its sample function is almost surely continuously 
differentiable on [a,b] of order k. 


( G(t — u, x — y) — G(—u , x — y)) W H (du, dy). 


Proof : The method of proof is similar to those of [T2j Proposition 3.1] and [411 Theorem 
4.8], but is more complicated in our fractional-colored noise case. 

The mean square derivative of Y at t € (0, oo) can be expressed as 

Y '(t)=[ [ G'(t-u,x-y)W H (du,dy), 

J — oo J 


where G' := dG/dt. This can be verified by checking the covariance functions. For every 
s,i £ ( 0 , oo) with s <t, similarly to the proof of Theorem [T] we derive 


E(|y'(t) - T'(s)| 2 ) = E (7° / ( G , (t-u,x-y)-G , (s-u,x-y))W H (du,dy)) 

\J — oo JJ 


oo poo 


\(\‘-^rr v ^ dv 

4(2v r) d J R d J o J o \u-v_ 

x ( e -\(t+nm 2 _g-|p+4l«| 2 ^ ( e -^t+v)\tf _ e -I( s+ „)|e| 2 


2-2 H 


In the above, we have used the fact that the Fourier transform of the function y i —> G'(t+u, y) 
is 

§- t (FG(t + u,-m) = . 

Denote by J-o,oo the restricted Fourier transform of / £ L 1 (0, oo) defined by (-Fo.oo/Xt) = 
J 0 °° e~ lXT f(x)dx, r£l, By applying the Parseval relation ([9]) for the restricted transform 
(see Lemma A1 in 0), we see that for all s,t € [a, b] with s < t, 


E(|y'(t) 
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Hence, as in nasi], by using Kolmogorov’s continuity theorem, we can find a modification 
of Y such that Y(t) is continuously differentiable on [a,b\. Iterating this argument yields 
the conclusion of Theorem [21 ■ 


Remark 2 In J23]/, the authors obtained a similar decomposition for the solution to the 
linear heat equation with white noise in time and Riesz covariance in space in dimension 
d = 1 . Our Theorem [7] shows that the time behavior of the process u is very similar to the 
behavior of the bifractional Brownian motion (see the main result in mb 


By applying Theorems UJ [2] and the results on uniform and local moduli of continuity 
for Gaussian processes (see e.g. [20l Chapter 7] or [22]), we derive the following regularity 
results on the solution process {u(t,x),t > 0}, when x E W l is fixed. For simplicity, we 
avoid the point t = 0. 


Proposition 1 Let be fixed. Then for any 0 < a < b < oo, we have 


lim 


sup 


| u(t, x) — u(s, x)| 


= C 3 a.s. 


e ^° S ,te[a,fc],| S -t|<£ |s - t|T'y / log(l/(t - s) 
where 0 < C3 < 00 is a constant that may depend on a, b , 7 and x. Or 


lim 


I u(t, x) — u(s, ®)| 
sup - :— = C3 a.s. 


£ ^-°s,tela,b],ls-tl<£ e'lyjlogil/e) 

Here and in Propositions [1 and [3J 7 = H — . 

Proposition 2 Let x E M. d be fixed. Then for any to > 0 we have 

— sup, t _ tot<£ | u(t, x) - u(t 0 , s)| 

lim- — -= ca a.s., 

£ ^° e 1 V lo g log(l/e) 

where 0 < C4 < 00 is a constant. 


By applying Theorems [TJ [2] and the Chung-type law of iterated logarithm in [23], 
we derive immediately 


Proposition 3 Let x E M fZ be fixed. Then for any to > 0 we have 


sup| f _i 0 |< £ \u(t, x) - u(to, x)\ 
TTO (e/loglog(l/e))7 


a.s., 


where 0 < C5 < 00 is a constant depending on the small ball probability of the Gaussian 
process U in Theorem\]\ 


Further properties on the local times and fractal behavior of the solution process 
{u{t,x),t > 0}, when x E is fixed, can be derived from [ 551 [39l HD] . 


10 









2.2 Sharp regularity in space 

In this section we fix t > 0 and analyze the space regularity of the solution {u(t, x),x € M. d }. 
We start with the following result. 


Theorem 3 For each t > 0, the Gaussian random field {u(t,x),x € M rf } is stationary with 
spectral measure 


A (d£) = a H {2-K 


\—d 


fj 


dudv 


«)ier 


o Jo \u-v 


2-2 H 


h(di). 


Proof : It follows from the Fourier transform of the Green kernel in m and Parseval’s 
identity Q that the covariance function of {u(t,x),x € M d } is 


E(u{t,x)u(t,y)) 



The conclusion of Theorem [3] follows. ■ 

It has been shown in [31 Proposition 4.3] that there exist two strictly positive con¬ 
stants c\,h,C 2 ,h such that 


Cl,H(t 2H A 1) 



< 


// 


dudv 


(f+f)ier 


\u — v 


2-2 H 


< C2,H(t 2H + 1 ) 


1 + ICI 2 


2 H 


(27) 


This and Condition m imply that the spectral measure A (dfi) is comparable with an 
absolutely continuous measure with a density function that is comparable to |£|-( Q + 4 ^) 
for all ( £ ^ with |£| > 1. As shown in [29l 19, [39], [40] . this information is very useful 
for studying regularity and other sample path properties of the Gaussian random field 
{u(t,x),x € M. d }. In the following we show some consequences. 

We start with the following estimate on E(|w(i,a;) — u(t,y) | 2 ). To this end, let 
ft = min{l, 2 H — and let 

p = { 1 n = 1 ' 

I 0 otherwise. 


Theorem 4 Assume that [T9 1 ) and V22\) hold. For any M > 0 and t > 0, there exist positive 
and finite constants cq, c-j such that for any x, y € [—M,M] d , 

c&\x - y\ w [log < IE( \u(t, x) — u(t, y)\ 2 ) < c 7 \x - y\ 2p [\og — l _ . (28) 
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Proof : Take i,i/6 [— M, M] d and let z := x — y E M rf . Using again (flTT) and Parseval’s 
identity, we can write 


E( | u(t, x + z) — u(t , x)\ 2 ) 

. . w /" f f* dudv 

= a H (27r)- d 



o Jo \u~v 


2-2 H 


.(l-cos <£,*»§; f f 


I 


e ' — 1 e 2 e 2 


dudv 


\u — v 


2-2 H 


Mr 

e 2 


Let us first prove the lower bound in (1251) . Using (12711 . with a generic strictly 
positive constant depending on t. H (that may change from line to line) and the lower bound 
in (1221) . we derive 


E | u(t, x + z) 


u(t,x) I 2 


- Cl,H '‘Xi>ilfF (i+wi (1 - cosg ’ 2>) 

- Cl ^L,lfF® (1 - cosfc2>) ' 


By making the change of variables using spherical coordinates, we have 


E |tt(f, x + z) 


ti(t,x)i 2 > ci j// ,tNr d+a+4 - f/ 



.d-i-a-4H( 1 _ cos (r(0,9 z ))dra(d0), 

(29) 


where = |f| and a(d6) is the uniform measure on the unit sphere E> d 1 . 

Next we distinguish three cases: (i) 2 H — > 1, (ii) 217 — = 1 and (iii) 

2H - 4=2 < 1 

In case (i) and (ii), we observe that for |z| small, 



,«J-l-a-4H( 1 _ cos ( r .^, 0a )) dr . 


> J I" r d - l -«~ m (r(9,9, z )) 2 dr 
2 -'Ll 

> C|z| d+ 2 -“- 4 // (log^) P ( 0 , 0 2 ) 2 , 


where the extra factor log ^ appears in case (ii). Plugging this into (1291) gives the desired 
lower bound. 

In case (iii), the integrand r 1 —> r d ~ l ~ a ~ iH (1 _ cos (r(9,9 z )) is integrable at both 0 
and infinity. The fact that 1 , 1/6 [— M, M] d ensures the integral has a positive lower bound. 
This gives the lower bound in (128j) . 

Now we verify the upper bound in (|28ji . Similarly to the above, the right-hand side 
of (1271) and (1221) imply 

E( \u(t, x + z) — u(t, x)\ 2 ) 

< c 2 ,H,t\z\- d+a+4H [ [°° r d ~ l ~ a ( , (1 - cos(r( 0 , 9 z ))dra(d9). 

Jsd-lJo \\A + r J 
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Again, by distinguish three cases: (i) 2H— > 1, (ii) 2= 1 and (iii) 2 < 1, 

we can verify that the upper bound in (12811 . Since this is elementary, we omit the details. 


Theorem U] suggests that the sample function x i—)• u(t,x ) is rough (or fractal) when 
2 H — < 1, and is differentiable when 2 H — > 1. This is indeed the case as shown 

by the following theorem. 


Theorem 5 Assume that [TV 11 and [22 1) hold and t > 0 is fixed. If 2 H — > 1, then 

{u(t,x),x € has a modification (still denoted by the same notation) such that almost 
surely the sample function x H > u(t,x ) is continuously differentiable on R“\ Moreover, for 
any M > 0, there exists a positive positive random variable K with all moments such that for 
every j = 1 ,... ,d, the partial derivative g§r u(t,x ) has the following modulus of continuity 
on [-M,M] d : 


sup 

x, y £ [—M , M ] d , | x—y \ <e 


d . . 


d , . 

-K—u{t,y) 

dyj 


< Ke 2H -^ 


-l 



(30) 


Proof : The method of proof is similar to that of Theorem [2] above or [411 Theorem 4.8]. 
By applying Theorem [31 we can show that the mean square partial derivative -u(t, x) 
exists and 


E 


d . 

dT, uit ' x) 


d 2 

-^—u(t,y) 

dyj 


f 

J R c 


A(d£) 

r J ' 

co S (f,x- y)) 


< C \x - y\* H+a ~ d - 2 . 


In the above, we have used the fact that 0 < 4 H + a — d— 2 <2, so the last inequality can be 
derived as in the case of fractional Brownian motion. Finally the desired result follows from 
the general results on modulus of continuity of Gaussian random fields (cf. [20, Chapter 7] 
or [50j Section 4]). ■ 

Finally, we consider the non-smooth case. For simplicity, we assume 2 H — < 1. 

The case 2 H — = 1 is more subtle and will require significant extra work. 

By combining Theorem [3] and relation (1271) with the results in [29] Section 8 ] (see 
also mm), we obtain the following useful lemma. 

Lemma 1 Suppose 2 H — < 1. Then ,for every fixed t > 0, the Gaussian field 

{u(t,x),x € M d } is strongly locally nondeterministic. Namely, for every M > 0, there 
exists a constant cs > 0 (depending on t and M) such that for every n > 1 and for every 

x,yi,-,y n e [-M, M] d , 

Vai(u(t,x)\u(t,yi),...,u(t,y n )) > eg mm {\x - yj\ 4H+a ~ d }, 

0<j<n 

where yo = 0 . 
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Because of this, the Gaussian random field {u(t,x),x E M d } shares many local 
properties with fractional Brownian motion B a = {B^(x), x E M d } with /3 = 2 H — We 
give some examples. 

Proposition 4 (Uniform and local moduli of continuity) Suppose 2 H — < 1. Let t > 0 

and M > 0 be fixed. Then 

(a) almost surely, 

max l6 d ,Ift|<e I u(t,x + h) - u(t,x)\ 

inn--- - = cq. 

sPy/\og(l/e) 


(b) For xq E M. d , 


max,M< £ | u{t, x 0 + h)~ 

lim- L -^—-— - 

£ -> 0 eh ydog log(l/e) 


u(t,x 0 )\ 


Go- 


In the above, /3 = 2 H — and 0 < cg,cio < oo are constants. 


Proof : The conclusion (a) follows from Lemma Q] and [221 Theorem 4.1], while (b) follows 
from [22] Theorem 5.1]. ■ 

Proposition 5 (Chung’s LIL) Suppose 2 H — < 1. Then for every t > 0 and xq E M. d , 

we have 

max| fe |< e \u(t,x 0 + h) - u(t,x 0 )\ _ 

£^/(log log(l/e)^ Cl0 ’ 

where (5 = 2 H — and 0 < cio < oo is a constant. 


Proof : It follows from Lemma |T] and Li and Shao HZ]. See also [J8] Theorem 3.2 ] for 
related results. ■ 

We conclude with the following remarks. 

Remark 3 (i) Propositions [^] and{5\ show that, when /3 = 2 H — < \ an d t > 0 is 

fixed, the solution process {u(t,x),x E M d } behaves locally like a fractional Brownian 
motion {B^{x),x E M d }. This is not surprising. Indeed, similarly to Theorem 0 and 
m Proposition 3.1], one can show that {u(t,x),x E M^} is a smooth perturbation of 
{Bh( x ),x E R d }. 

(ii) Further properties on the local times and fractal behavior of the solution process 
{ u(t , x),x E M d } ; when t > 0 is fixed, can be derived from T3R I ~39[ It is also pos¬ 

sible to investigate sample path properties of the Gaussian random field {u(t,x),t > 
0, x € M rf } in both time and space variables. The methods developed in m will be 
useful for this purpose. We will study these problems in a subsequent paper. 
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